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ABSTRACT 



3/008 


An approximate theory is developed to find out whether a star, moving initially on the plane of sym- 
metry of a galaxy, escapes during the collision with another galaxy of small dimensions. The escaping 
stars are inside an ellipse on the plane of symmetry of the galaxy. Numerical calculations in a Bottlinger 
model of a galaxy give results of the same order of magnitude. The total amount of the escaping matter 
is only a few per cent of the total mass of the galaxy, if the collision velocity is of the order of 2000 
km/sec. The behavior of the escaping and non-escaping stars is described. The escaping stars do 
follow the colliding galaxy and do not form a bridge or link between the two galaxies. 

I. ESCAPING STARS 


2000 
o not 


The purpose of this paper is to consider the effects of the collision of two galaxies on 
stars moving initially on the plane of symmetry of one galaxy. 

If we assume a relative velocity of collision of about 2000 km/sec (Spitzer and 
Baade 1951), the collision lasts about 10 7 years. During that time the perturbations due 
to the ellipticity of the first galaxy are insignificant for stellar orbits close to its plane of 
symmetry (Contopoulos and Bozis 1962; hereafter referred to as “Paper I”)- Therefore 
in order to study approximately the effects of the collision on these orbits we may use a 
spherical model of a galaxy. The colliding galaxy is also assumed to be spherical but of 
small dimensions. We will see that, if its diameter is less than 2 kpc, its effects are es- 
sentially the same as those of a point mass. In the numerical examples below the mass of 
both galaxies is assumed to be about 10 u M o, and the radius of the first galaxy is taken 
10 kpc. 

In Paper I it has been proved that, if two such galaxies collide with a relative velocity 
v = 2000 km/ sec, their relative motion may be considered as uniform; in fact the change 
of velocity is of the order of 30 km/sec and the change of the direction of the velocity 
is less than 5°. 

Further, during the collision, a star moving in a circular orbit on the plane of sym- 
metry of the first galaxy 7 O with radius 8 kpc describes only about 1/20 of a revolution, 
i.e., about 18°; therefore, we may 7 consider, in a rough approximation, the unperturbed 
orbit of A to be a straight line, x' = x 0 ', perpendicular to the axis x', which is directed 
away from the center O (Fig. 1), and parallel to the axis y ' . Under these assumptions the 
perturbation may be considered as a two-body encounter that lasts as long as the 
collision itself. 

The galaxy O' is assumed to move along the axis s' with uniform velocity 7 v , reaching 
the point Oi on the plane xOy at t = t 0 . The initial velocity of the star is v c ; if the per- 
turbing galaxy was not present, S would reach the axis x' at t — t\. 

The velocities v and v c are of the order of 2000 km/sec and 200 km/sec, respectively; 
therefore, during 10 7 years the galaxy O' moves approximately from s' — 10 kpc (initial 
position) to s' = — 10 kpc (final position), whiled moves about 2 kpc along a parallel to the 
y'-axis. The influence of O' on S when the distance OO' is larger than 10 kpc is very 7 small. 
In fact, as we will see below, the escaping stars are inside an ellipse whose dimensions are 
of the order of 4 kpc. Therefore, their minimum distance from O' is about 2 kpc or less. 
When a point comes from infinity to a distance of 10 kpc the relative velocity 7 changes 
only 7 by about 1 per cent and the change of its direction is of the order of 2 per cent of 
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the angle between the asymptotes. In fact, if v m is the relative velocity before the en- 
counter and v the relative velocity at a distance r from O', then 


hence 


9 2 Gm , 

v 2 — \~ v 


Av = v — v 


Gm 


and for r = 10 kpc, = 2000 km/sec, and m = 10 u Mo, 

Av = 2 2 km/sec. 



Fig. 1. — A galaxy O' of mass m is moving downward along the axis 2' toward the point Oi on the plane 
of symmetry of the galaxy O, of mass M. The axis Ox is parallel to the projection 0\X n of 0\z' on the 
plane of symmetry. The distances of O' from O, Oi, and S are d* } d, and D , respectively. 


The angle 2p between the initial and final velocities is given by tan p — p/q, where 
Xi ~/p 2 — yi 2 /g 2 = 1 is the equation of the relative orbit. As q is of the order of 2 kpc, 
and p = Gm/ vj 2 0.1 kpc, we have p/q ~ 0.05. 

The tangent at a point {x\, yi) whose distance from the focus is 10 kpc makes an 
angle p r with the axis yi, where 

tan^'=tan^(l-^), 

and [*! - V ( p 2 + <? 2 )] 2 + yr = 100. 

Then the deviation (i p — ip') is given by 


. (p/q)[10 + p- V(100 + 20j>- g 2 )] 

^ V) 10 + p+(pi / q2) V ( 100+20/.- 9 2 ) 
Hence 

\f/ — 1 p r 

^-^-=^0.005 q 2 , 

1 




p 0.05 q 2 
n' 10 


= 0.005/>g . 
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and for q = 2, this is 2 per cent. Therefore we may consider the initial and final relative 
velocities of S with respect to O as being equal to the corresponding velocities in a two- 
body encounter at infinite distances. 

The initial relative velocity of the star S with respect to the galaxy O' is t’o — v c — v; 
the corresponding final velocity Vq is equal to vq in measure but not in direction. The 
final velocity of S with respect to the galaxy 0 is v/ = v Q f + v; if v f is greater than the 
velocity of escape v e , the star escapes from the galaxy O. 

In what follows we use a frame of reference where O' is at rest. The galaxy 0 is moving 
with respect to O' with velocity — v along the axis z' that makes an angle e with the 
plane xOy (Fig. 1). The projection of z' on the plane xOy makes an angle 9 with Ory', 


and <p is the angle z'Oiy'. Then 

cos ^ = cos e cos 9 . 

( 1 ) 

At the initial time t we have 

y' = OiS' = vc(t-h) 

( 2 ) 

and 

z' = OiO'= v(k— /). 

( 3 ) 

If a? is the angle between t> 0 and — v, we have 


sin co __ sin tp _ sin ( <p — co ) 

( 4 ) 

v e vq v 

and 

Vq 2 = V 2 + Vc 2 +2 vv c cos <p . 

( 5 ) 


The corresponding angle between v 0 ' and — v is co', and the angle between and vj 
is 2\f/, which is approximately equal to the angle between the asymptotes of the hyper- 
bolic orbit. Therefore 


tan = 


Gm 

qv o 2 ’ 


( 6 ) 


where q is the impact parameter. 

The impact parameter is found as follows. Let SS' be parallel to the axis x ' and S'Q' 
parallel to v 0 , O'Q and O'Q' perpendicular to SQ and S'Q', respectively; S'Q' is on the 
plane S'OiO' and QQ' is perpendicular to SQ and S'Q'. The angle between the planes 
QQ'S' and O'Q'S ' is n, and the angle between the planes QO'S and Q'O'S' is v. The angle 
Q'S'Oi is equal to <p — u, therefore 

q' ~ — y' sin (<p — co ) — s ' sin co = — ^y'+z' sin ( tp — to ) 


= — *o ) sin ( (p — a; ) = — y</ sin ( <p — a? ) , 


( 7 ) 


where 


V = v c (1 Q —li) 


( 8 ) 


is the y' coordinate of S at / — to. The term q' is positive if yo < 0, and negative if 
yo > 0. Further 

S'S — Gi0 2 = xo', Q , Q=x 0 / sin5, and q 2 = q' 2 + x 0 ' 2 sin 2 5 — Iq'xo' sin5 cos n . ( 9 ( 

From the trihedral angle S'Q'SOi we find 

cos 5 — sin ((p — co ) cos I (io) 

sin S cos n= — cos ( <p — co ) cos I 


and 


(it) 
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where 7 is the angle between the planes y'Oiz' and y'Oix', hence 

cos e sin0 


cos 7 = cot <p tan 0 — 


Therefore, 


sm <p 
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(12) 


q 2 = x Q ' 2 [ 1 — sin 2 ( <p — co ) cos 2 I ] + y 0 ' 2 sin 2 ( <p — co ) 
— 2xo'yo ' sin (<p — co)cos(<p — co)cos 7 . 

A star escapes if 

V/ 2 — v 2 -\-Vq 2 —2vv 0 cos c o' > z> e 2 . 

This is possible only if v + v 0 > »«. Then we can set 


-1 < 


^ 2 +^ 0 2 — V e 2 ,^V 2 +Vq 2 —V c 2 V+V c COS<p ^ . 

= cos co 0 < = = cos CO < 1 , 


( 13 ) 


( 14 ) 


( 15 ) 


( 16 ) 


( 17 ) 


2 Wo ^ wu ' 2 Wo vo 

because v c < v e . Therefore 0 < co 0 ' < t, and if a star escapes we have 

cos co' < cos co 0 ' . 

The angle co' can be found from the trihedral angle SjVoVq'(—v); 

cos co' = cos 2\p cos co + sin 2 ^ sin co cos ( tt — i>) , 

where (ir — v) is the angle between the plane S/VqV 0 ' (i.e., the plane O'QS of the relative 
orbit) and the plane S/Vq(— v) (which is parallel to the plane O'Q'S'); thus (ir — v) is 
the complement of the angle QO'Q ' and 

q' — xo sin 5 cos ju , . . , , , 

cosy = = i [ — y 0 sm ( c? — co) + Xo cos(<p — co) cos 7] . (18) 

q 9 

If we use relations (6) and (17) in inequality (16) we get 

{ Xq' 2 [ 1 — sin 2 (<? — co)cos 2 7] + y 0 ' 2 sin 2 (^ — co) — 2%o Jo sin(<p — co)cos(<p — co)cos7 } 

, , 2 Gm . . . . . 2 Gm . , N r , 

(cos co — cos co 0 ) H y sm 00 sin (^? — co) y 0 y Sin w cos {<p— co) cos Ixq 09) 

V (p 2’ o 

G 2 m 2 

< Y ( COS CO + cos COq / ) , 

i’o 4 

i.e., the escaping stars are inside the ellipse 

Xo' 2 ( 1 — sin 2 ( <p — co ) cos 2 7 ) + ( y 0 f — y c ' ) 2 sin 2 ( ^ — co ) 

— 2xq'( yo — y c ')sin (<p— co ) cos (<p— co ) cos 7= — ^ m S * n — ~ — — 

' ' ^ • fl 0 4 (cosco — COSCOo ) 

whose center is {xo — 0, y c '), where 

, _ Gm sin co _ 2 Gmv c 

Vo 2 si n( v? — co)(cos co — cos co 0 ') ( v e 2 — fl c 2 ) t'o* 

The axes of this ellipse form an angle y with the axes x' and y\ where 

— tan ( <p — co ) cos e sin 0 


M2’ 


( 20 ) 


( 21 ) 


tan y ~ — tan ( <p — co ) cos 7 — - 


sin 


( 22 ) 



ESCAPE OF STARS 


1243 


No. 4, 1964 


For 0 < 0 < x, the angle y is positive (measured counterclockwise) if <p — co > 7r/2 
and negative if ^ — co < x/2. For — tt < 6 < 0, i.e., if the projection of s' on the plane 
x r O\y f is on the negative side of the axis x\ the angle y is positive if <p — co < tt/2 and 
negative if ip — co > tt/2. For 6 — 0 or 8 = ir, the angle 7 = 0. 

If 8 = 7r/2 then = tt/2, / = e, and equations (20)-(22) become 

x 0 ' 2 ( 1 —cos 2 co cos 2 e) + cos 2 co(y 0 ' — y/) 2 — 2 jc 0 ' ( y 0 f — y c ) sin co cos co cose 

= G ~™ 2 sin2 “V (23> 

t'o 4 (cos £0 — cos coo') 2 ’ 


and 


Gm tan co 

I’o 2 (cos co — COSCOq') 1 


tan 7 = — cot co cos e . 


( 24 ) 

( 25 ) 


From equation (20) we find easily the area of the ellipse 


F _ , irG 2 m 2 sin 2 co 0 ' 

7ra * a v i’o 4 (cos co — coscoq ') 2 sin ( ^> — co ) sin I' 


( 26 ) 


TABLE 1 


r 

(kpc) 

a x | 

(kpc) 

ay 

(kpc) 

E 

(kpc) 2 

yc 

(kpc) 

o) 

(degrees) 

4 

1.41 

1.42 

6 

—0.7 

6.22 

6 

1.94 

1.95 

12 

-1.2 

6.55 

8 

2.41 

2.42 

18 

-1.6 

6.19 

10 

2.80 

2.81 

25 

-1.9 

5.74 


When t>— ► oo 5 then E — » 0. If I — i r/2, then the ellipse becomes 

,0 , , , ^ G 2 w 2 sin 2 co 0 ' 

Xq 2 + ( y 0 — y c ) 2 sm 2 ( (p — co ) = — yrz . 

w ^O 4 (cosco — COSCOo ) 2 

The semi-axes of the ellipse are then 


a x = 


Gm sin co 0 ' 2Gm sin co 0 ' sin ( <p — co ) 


and 


r 0 2 (cos co — cos coq') (i' e 2 — v c 2 )s[n <p 


Gm sin co 0 ' 


2Gm sin co 0 ' 

^o 2 (cos co — cos co 0 ')sin (<p — co) ( v e 2 — fl c 2 )sin <p 


If v — > 00 , then — > 0 and a* — > 0. 

The ratio of the axes is 

a x /a y — sin(<p — co) = v sin <p/ v Q , 
and the area of the ellipse is 


E = ira x a y = 


47 rG 2 m 2 sin 2 co 0 ' sin ( — co ) 
( Ve 2 — ?c 2 ) 2 sin 2 <p 


( 27 ) 


( 28 ) 


( 29 ) 


( 30 ) 


( 31 ) 


Table 1 gives the values of a x , a y , 2£, y c ', and co for the model considered in Section II 
with e = tt/2, <p = tt/2, = 2000 km/sec = 2.0454 kpc/10 6 years, and different mini- 

mum distances r ~ V(a 2 + £ 2 ) of the colliding galaxy from the center. 
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If the effective radius of the galaxy R G is about 10 kpc, the above area is between 2 
per cent (for r = 4 kpc) and 8 per cent (for r — 10 kpc) of the total area. 

The proportion of the escaping stars is roughly Ep/wRcp 2 where p is the mean density 
of the galaxy and p is the mean density in the region of escape. 

Table 2 gives the values of <p, E , y c \ co for (e = 75°, 6 = 0°), (e = 45°, 6 = 0°) and 
(e = 45°, 6 = 45°). 

The value of y in the first two cases is zero, and in the last case varies between 
— 39?8 and — 39? 1 . 

It is seen that the values of E and y/ in all three cases are quite near the corresponding 
values of Table 1. Therefore, the amount of escaping matter from the plane of symmetry 
of a galaxy does not depend very much on the direction of the collision. 

We do not give the values of E , etc., for r < 4 because near the center the assumption 
of a binary collision between the star and the approaching galaxy is not even approxi- 
mately satisfied. 

If the diameter of the colliding galaxy is 2 kpc (or less), its cross-section is small, in 
general, in comparison with the area E of the region of escape; therefore, the amount 
of escaping matter is essentially the same as in the case of a point mass. If the diameter 
is larger than 2 kpc, the amount of escaping matter will be smaller. 


TABLE 2 



ii 

•-j 

(SI 

(t> 

II 

O 

= 75° 

e = 45° t 0 = 0°, <p 

= 45° 

e = 45 

e=4S°, * 

> = 60° 

r 

(kpc) 

£ 

(kpc) 2 

y c 

(kpc) 

CO 

(degrees) 

E 

(kpc) 2 

yc 

(kpc) 

to 

(degrees) 

£ 

(kpc) 2 

y c 

(kpc) 

(degrees) 

4 

6 1 

-0.7 

5.85 

7 ! 

-0.7 ! 

4.10 1 

7 1 

-0.7 

5.12 

6 

11 ! 

-1.2 

6.15 i 

11 

-1.1 

4.29 ! 

13 

-1.2 

5.37 

8 

17 

-1.6 

5.82 

17 

-1.5 

4.08 

20 

-1.5 

5.09 

10 

23 

-1.9 

5.41 

23 

-1.8 

3.80 

28 

-1.9 

4.74 








II. NUMERICAL CALCULATIONS 


In order to calculate numerically the orbits of stars during the collision of two galaxies 
we have chosen a simple model for the first galaxy, in which the forces per unit mass are 
given by Bottlinger’s fonnula 


a 


1 + br*' 


(32) 


This formula has been applied by Lohmann (1953; 1954a, b) to our Galaxy, as well as 
to M31 and M33, and by Kerr" and de Vaucouleurs (1956) to the Large Magellanic 
Cloud. It is discussed, together with other models, by de Vaucouleurs (1959) and Perek 
(1962). 

The values of a and b were chosen to represent approximately the force function F 
given by Schmidt (1956) for the plane of symmetry of our Galaxy, namely, 

a= 0.0055 (10 6 years)- 2 , b= 0.012 kpc" 3 . 


The corresponding potential function is 


U = 



r 2 — kr + k 2 l 

(r + k)* J 



where k — Zu 1/3 and U « 


( 33 ) 


= 0 . 
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In this model the mass 

of the Galaxy is 


The circular velocity is 

M = 1.02 X 10 l W o . 

bG 

( 34 ) 

and the escape velocity is 

I'c = V ( Fr ) , 

( 35 ) 

t'e = V(2U). 

( 36 ) 


Table 3 gives the values of v c and v € for different distances r from the center in the 
above field. 

Now we consider the following situation (Fig. 1): A point mass O r {m) goes through a 
spherical galaxy 0(M), represented by the force function (32). This force function is as- 
sumed to be stationary; this is justified if the colliding galaxy does not come very close 

TABLE 3 


r 

(kpc) 

Vc 

(kpc/10 6 years) 

Ve 

(kpc/10 6 years) 

f ! 

(kpc) 

(kpc/ 10 s years) 

Ve 

(kpc/10 6 years) 

4 

0.2231 

0.4294 

14 

0.1783 

0.2549 

6 

.2348 

.3750 

16 

.1676 

.2388 

8 

.2220 

.3321 

18 

.1584 

.2253 

10 

.2057 

.2997 

20 

0.1506 

0.2138 

12 

0.1909 

0.2748 





to the center of the galaxy O, because then the amount of escaping matter is rather 
small. In fact the density corresponding to the force function (32) is 

3a 

P ~ 47rG(l + 6r 3 ) 2 ’ 

and the mean density is 

- _ 3a 

p ~THGbR?' 


so that if we take a mean radius of the galaxy R G — 10 kpc, then the proportion of the 
escaping stars is 

Ep 0.04E 

7rR G 2 p~ (l + 0.012r 3 ) 2 ’ 


for r = 10 kpc this is about 0.6 per cent while for r — 4 kpc this is 8 per cent. 

The motion of O' is assumed to be uniform along the line O f O \ , where Oi(a, fi) lies on 
the plane xOy. The plane xOz is taken parallel to the direction O'Oi and e is the angle 
between this line and the plane xOy. The coordinates of the star S are x, y, s, and d , 5*, 
D are the distances of O' from Oi, O, and S, respectively. 

The unperturbed orbit of 5 is, in general, a rosette on the plane xOy . In this paper, 
however, we consider only the case of circular orbits around 0. 

The perturbation acting on the star is equal to the acceleration of S due to the galaxy 
O' minus the acceleration of the galaxy O, taken as a whole. As the attraction between 
O and O' is mF(d *), the acceleration of 0 is ( m/M)F(d *). Thus the equations of motion 
of S are 


d 2 x 

If 

d 2 y 
dt 2 


= ~F(r*)— *+^r{a + d cose- x) -F{d*) 


D 

Gm 


m (a + d cos e ) 
M d * 


( 37 ) 


m P 




M d* 


(38) 
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and 





d 2 z __ 
dt 2 * 


d 

sin e 
d* 1 

( 39 ) 

where 






r*= V(x 2 +y 2 + z 2 ), 



( 40 ) 


d*= V [ ( a + cos e) 2 +/3 2 + (d sin e ) 2 ] , 



( 41 ) 


D = V [ ( a + d cos e — * ) 2 + ( j3 — y ) 2 + ( 

d 

sin e — z ) 2 ] , 

( 42 ) 

and 




Pi- 

ll 

o 

1 



( 43 ) 


The numerical values used are the following: m/M = 1, MG — a/b , v = 2000 km/ 
sec = 2.0454 kpc/10 6 years, do = 100 kpc, a = 8 kpc, /3 = 0, and e = 90° (case A), or 




Fig. 2. —The position (r, g in polar coordinates) of escaping ( open circles) and non-escaping {solid 
circles) stars in their unperturbed orbits when the colliding galaxy is going through the point 0/(a = 
8 kpc, /3 = 0). In case A, e = 90°; in case B, e = 45°. 

e = 45° (case B). The initial conditions for i = 0 were taken in such a way that the star, 
moving in the unperturbed circular orbit, would reach an axis forming an angle g with 
the positive x-axis (g is positive if measured counterclockwise), when the colliding galaxy 
reaches the point Oi (a, /3), namely, at about t = 48.89 X 10 6 years. 

The greatest part of the calculations was made with the IBM 7090 computer of NASA's 
Institute for Space Studies, New York. The Runge-Kutta integration method was used 
with a step 0.1 (10 6 years) and the orbits were calculated for 200 X 10 6 years. At t = 
200 X 10 6 years the colliding galaxy is at a distance of about 300 kpc from the galaxy O 
and its influence is then negligible. In some cases, when the stars escaped rather fast from 
the galaxy O, the calculations were stopped earlier. 

The initial z coordinates and velocities were all assumed to be zero. Figure 2.4 
shows the positions that the stars would have in case A (e = 90°) at t = 48.89 X 10 7 
years (when the colliding galaxy reaches the point Oi), if there was no perturbation. The 
polar coordinates of these positions are (r, g); therefore the initial polar coordinates, for 
1 = 0, are 

r 0 = r 


go = g ~ 48.89 — X 5 7?295 7 795 

r 


and 
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The dots represent stars that are not escaping. The small circles represent escaping 
stars. The region of escaping stars is approximately a circle of diameter 5 kpc. From 
Table 1 we see that the ellipse of escaping stars for r — y/{a 2 + 0 2 ) = 8 kpc is very nearly 
a circle of diameter 4.8 kpc; the agreement with the calculations is very good in this case. 

The center of the circle in Figure 2.4 is not at x c ' = 0, as it was expected theoreti- 
cally, but at about x c f — 1 kpc, i.e., further out than the point Oi. On the other hand, the 
y f coordinate of the center is about yj = — 1.5 kpc, i.e., very close to the calculated value. 
Therefore, the only discrepancy between theory and calculations is that the actual region 
of escape is further out than expected; this is due, of course, to the fact that the unper- 
turbed orbits are curved and not straight lines as they were assumed in Section I. The 
only consequence of this effect is that the region of escape actually includes less dense 
areas, therefore, the amount of escaping matter is smaller. 

In case B (e = 45°) the calculated orbits are marked in Figure 2 B. The region of 
escape is approximately an ellipse of dimensions 2a — 8 kpc, 2 b ~ 5 kpc, and the area 
E — 31 kpc 2 . This is bigger than expected on the basis of the simple theory developed in 
Section I. However, even in this case the area E is only 10 per cent of the effective area 
in the plane of symmetry of the galaxy 0, and as it includes mainly the outer parts of 
the galaxy the mass of the escaping matter is less than 2 per cent of the total matter 
near the plane xOy. 

The individual orbits of stars during the collision show a number of interesting char- 
acteristics. Figures 3 and 4, copied from plots made by the NASA computer, show the 
main t\pes of orbits encountered, in the case e = 90°. 

A common characteristic of all the orbits is that the perturbations are very small until 
the colliding galaxy goes through the galaxy O. Until that moment the orbits are prac- 
tically circles on the .yy-plane. After the collision the orbits may be divided into three 
types: (a) orbits of stars that remain far from the region of escape; ( b ) “escaping” orbits; 
and ( c ) orbits of stars near the region of escape. 

If the star is far from the region of escape, its orbit is not changed much by the col- 
lision. It remains almost circular, and the s-motion is very small. Such is the case of 
orbit I (r — 8 kpc, g = —60°). 

The orbits of stars inside or near the region of escape can be divided into the following 
groups (cf. Figs. 2.4 and 3A-3D): (i) stars with g < 0° and r > 8, or r smaller than 8, 
but g not near 0°; (ii) stars with r < 8 and g < 0° but g near 0°; and (iii) stars with 
g > 0°. We discuss these groups separately. 

i) The stars with g < 0° that are inside or near the region of escape are accelerated 
toward the positive y'-axis as the colliding galaxy goes through the galaxy O. If r > 8 
kpc the orbits are deflected to the left. This deflection is bigger whenever the approach 
to the colliding galaxy is closer. For example, in case II (r = 9, g = 0°) the star is 
ejected abruptly to the left and down, while in case III (r = 11, g = —5°) the star 
escapes but with smaller velocity. In case IV (r = 11, g = 0°) the star does not escape, 
but its orbit does not differ much from the escaping orbits for t — 200 X 10 6 years. Only 
the fact that its velocity is less than the escape velocity makes it sure that it does not 
escape, but forms a very elongated orbit around the center of the galaxy O. 

All orbits near the boundary of the region of escape from (r = 11, g = 0°) to (r = 8, 
g — —25°) are similar to cases III and IV. The y component of the final velocity is 
negative for large r, while it is positive for r near 8. 

The same general behavior is shown by the orbits with r < 8 but g not close to 0°. 
For example, case V (r — 7, g = —20°) represents an escaping orbit; the star escapes 
toward the positive y direction and downward. 

ii) The orbits with r < 8 and g near 0° show an abrupt deviation to the right. This 
is most marked in case VI (r = 7, g = 0°) when the star is ejected abruptly to the right 
and downward. 

iii) The orbits with positive g are usually not escaping. In general, the additional 
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attraction of the colliding galaxy makes the orbits more elongated. Such is the case VII 
( r = 8, g = +30°), that is rather far from the region of escape. In case VIII (r = 8, 
g = +10°), however, the deviation is so large, that the orbit reverses the sense of rota- 
tion, although the star does not escape. A large deviation, but to the right, is shown in 
case IX (r = 7, g = +5°). Finally, case X (r = 8, g = +5°) shows an abrupt reversal 
of the motion of the star and ejection toward the negative y-axis. Such a behavior, how- 
ever, is quite exceptional. 



Fig. 3.™ Characteristic cases of orbits during a collision through Oi(a — 8 kpc, /3 = 0) with e = 90°: 
I (r = 8 kpc, g — —60°; non-escaping); II (r = 9 kpc, g = 0°; escaping); III (r — 11 kpc, g = —5°; 
escaping); IV (r ~ 11 kpc, g = 0°; non-escaping); V (r = 7 kpc, g = —20°; escaping); VI (r — 7 kpc, 
g = 0°; escaping); VII (r — 8 kpc, g = +30°; non-escaping) ; VIII (r — 8 kpc, g — +10°; non-escaping); 
IX (r = 7 kpc, g = +5°; non -escaping) ; X {r = 8 kpc, g = +5°; escaping). A, B give the projections of 
the orbits on the x^-plane, while C , D give the projections of some orbits on the as-plane. 


From the above discussion we conclude that the greatest proportion of the ejected 
stars escapes toward the left (between the positive y-axis and the negative x-axis) and 
downward. The stars do not follow the colliding galaxy after the collision and they do 
not form a bridge or link between the two galaxies. 

In case B (e = 45°) the orbits are similar to those of case A, except as regards the 
z coordinate (cf. Figs. 2 B and 4). 

i) The stars with r > 8, g < 0° have comparatively large motions upward, as the 
colliding galaxy approaches, which are changed to downward motions later, as in case 
XI (r = 14, g = 0°). This is because in these cases the approaching galaxy comes rather 
close to these stars before it reaches the plane xOy at Oi. 

If r is near 8 kpc but r > 8, and g ^ 0° the upward motions are so fast that the stars 
eventually escape upward and to the left (case XII, r = 9, g = 0°). If r is near 8 kpc 
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but g not near 0° the upward motions are small (case XIII, r = 8, g = —25°), and 
if r < 8 kpc the stars go downward, showing only a small upward kink, as 
in case V. The downward motions in this case can be explained easily, because as the 
approaching galaxy comes from the right, it does not influence appreciably the stars to 
the left of 0\ until it goes beyond the point Oi and comes rather close to them. 

ii) The stars with r < 8 but close to 8 and g near 0° show an abrupt motion down- 
ward (case XIV, r = 7, g — 0°). 



Fig. 4. — Characteristic cases of orbits during a collision through Oi(a = 8 kpc, 0 = 0) with e — 45°: 
XI (r - 14 kpc, g = 0°; escaping); XII (r = 9 kpc, g = 0°; escaping); XIII (r = 8 kpc, g = -25°; 
escaping); XIV (r = 7 kpc, g = 0°; escaping); XV (r = 9 kpc, g = +10°; non-escaping). Only the pro- 
jections on the x*js-plane are given. 


iii) The orbits of stars with g > 0° show all the peculiarities mentioned in case A. 
The z motions are of the form of case XV (r = 9, g — +10°) if r > 8, and of the form 
of case XIV if r < 8, but with less steep descent. 

We conclude that if a galaxy goes through the plane of symmetry of another galaxy 
at an angle different from 90° it causes the stars to acquire appreciable z velocities. The 
escaping stars move usually downward, but they do not follow the colliding galaxy in its 
course. 

The final conclusion of this paper is that, during a collision of two galaxies with 
rather large relative velocity, only a small fraction of matter escapes and this does not 
form any bridge or link between the galaxies. Therefore the observed bridges between 
pairs of galaxies probably are not due to collisions with relative velocities of the order 
of 2000 km/sec. It is possible, however, that much slower collisions, not considered in 
this paper, may cause much larger changes in the colliding galaxies and eventually ex- 
plain the formation of bridges between them. 
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